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A  computational  model  has  been  developed  in  order  to  simulate  the  thermal  and  electric  behavior  of 
thermoelectric  generators.  This  model  solves  the  nonlinear  system  of  equations  of  the  thermoelectric 
and  heat  transfer  equations.  The  inputs  of  the  program  are  the  thermoelectric  parameters  as  a  function 
of  temperature  and  the  boundary  conditions,  (room  temperature  and  residual  heat  flux).  The  outputs 
are  the  temperature  values  of  all  the  elements  forming  the  thermoelectric  generator,  (performance,  elec¬ 
tric  power,  voltage  and  electric  current  generated).  The  model  solves  the  equation  system  using  the  finite 
difference  method  and  semi-empirical  expressions  for  the  convection  coefficients. 

A  thermoelectric  electric  power  generation  test  bench  has  been  built  in  order  to  validate  and  determine 
the  accuracy  of  the  computational  model,  which  maximum  error  is  lower  than  5%. 

The  objective  of  this  study  is  to  create  a  design  tool  that  allows  us  to  solve  the  system  of  equations 
involved  in  the  electric  generation  process  without  needing  to  impose  boundary  conditions  that  are 
not  known  in  the  design  phase,  such  as  the  temperature  of  the  Peltier  modules. 

With  the  computational  model,  we  study  the  influence  of  the  heat  flux  supplied  as  well  as  the  room 
temperature  on  the  electric  power  generated. 

©  2009  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

We  have  developed  a  computational  model  that  simulates  the 
thermal  and  electric  behavior  of  a  thermoelectric  generator  that 
presents  some  advantages  with  respect  to  previous  ones. 

Some  applications  in  thermoelectric  generation  use  waste  heat 
fluxes  as  energy  source.  Benson  et  al.  [1]  investigated  geothermal 
heat  fluxes  with  temperature  gaps  of  80-180  °C;  other  gaps  are 
produced  in  the  ocean’s  currents,  in  solar  panels,  or  in  power  gen¬ 
erator  stations.  Matsuura  and  Rowe  [2]  propose  these  and  other 
residual  thermal  energy  sources. 

In  order  to  study  these  applications,  it  is  necessary  to  determine 
the  behavior  of  the  thermoelectric  devices.  The  analytic  solution 
most  used  to  solve  the  equation  system  of  the  thermoelectric  device 
is  the  ideal  thermoelectric  couple,  which  expression  was  obtained 
by  Ioffe  in  [3].  Other  references  that  use  this  expression  are  [4,5], 
and  in  the  case  of  thermoelectric  generation,  it  was  used  in  [6- 
1 1  ].  In  the  expression  of  the  ideal  thermoelectric  couple,  the  ther¬ 
moelectric  parameters  (Seebeck  coefficient,  electrical  resistivity 
and  thermal  conductivity)  are  not  dependent  on  the  temperature. 

A  better  approximation  is  to  consider  an  average  value  of  the 
Seebeck  coefficient  within  the  operating  temperatures  of  the  ther¬ 
moelectric  couple,  as  can  be  seen  in  the  researches  of  Ioffe  in  [3], 
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Arenas,  in  [12]  and  Yu  in  [13].  Nevertheless,  this  approximation 
does  not  adjust  to  reality  because,  in  fact,  the  Seebeck  coefficient 
varies  with  temperature. 

The  Buist  works  [14,15]  divide  a  thermoelement  into  segments 
and  apply  the  ideal  thermoelectric  couple  in  each  segment.  This 
model  keeps  in  mind  the  variation  of  the  thermoelectric  properties 
with  temperature.  Its  major  disadvantage  is  that  the  model  needs 
the  heat  flux,  the  hot  source  temperature  and  the  electric  current 
generated  as  input  parameters,  and  it  is  necessary  to  obtain  them 
experimentally. 

To  determine  completely  the  thermoelectric  device  behavior,  it 
is  also  necessary  to  consider  the  heat  exchangers  on  the  sides  of 
the  Peltier  module,  as  is  described  by  Stockholm  in  [16].  The  disad¬ 
vantage  of  Stockholm’s  model  is  that  it  simplifies  the  Peltier  mod¬ 
ule,  and  it  is  impossible  to  optimize  the  geometry  and  dimensions. 
Similar,  but  more  comprehensive,  models  have  been  posed  by  Are¬ 
nas  in  [12]  and  by  Kondraitev  in  [17],  but  these  models  do  not  cal¬ 
culate  the  transitory  state. 

Our  model  is  based  on  a  computational  model  developed  to 
simulate  the  behavior  of  thermoelectric  refrigerators  by  Vian  and 
Astrain  [18-20]. 

Our  model  solves  the  behavior  of  the  thermoelectric  generator 
and  makes  improvements  to  the  previous  models.  It  takes  into  ac¬ 
count  the  properties  of  the  thermoelectric  materials  as  functions  of 
temperature.  We  divide  the  Peltier  module  into  ten  different  parts, 
which  allows  us  to  determine  the  variation  of  the  thermoelectric 
properties  in  its  interior.  The  Thomson  effect  is  not  considered 
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Nomenclature 

A 

area  (m2) 

AT 

difference  of  temperatures  (K) 

C 

thermal  capacity  (J  K1) 

8 

test 

c 

specific  heat  (J  kg-1  K1) 

P 

density  (kg  m-3) 

Et 

electromotive  force  (V) 

Po 

electric  resistance  of  Peltier  module  (Q  m) 

I 

electric  current  (A) 

a 

thomson  coefficient  (VIC1) 

k 

thermal  conductivity  (W  m-1  K-1) 

T 

time  (s) 

L 

characteristic  length  (m) 

m 

load  ratio  ( RJRo ) 

Superscripts ,  subscripts 

N 

number  of  thermocouples  of  Peltier  modules 

c 

cold  side  of  Peltier  module 

P 

generated  power  (W) 

cont 

semiconductor  and  copper  contact 

a 

calorific  power  (W) 

Dissip 

dissipater 

Q 

calorific  power  per  unit  volume  (W  m-3) 

h 

hot  side  of  Peltier  module 

R 

thermal  resistance  (KW_1) 

i 

node  i 

Rcont 

contact  electric  resistance  (Q) 

in 

system  inlet 

R0 

electric  resistance  of  Peltier  module  (Q) 

k 

number  of  8  repetitions 

Rl 

load  electric  resistance  (Q) 

J 

joule  effect 

T 

temperature  (°C) 

j 

node  j 

T 

temperature  (K) 

leaks 

system  loses 

V 

temperature  in  next  time  step  (K) 

max 

maximum  value  of  studied  curve 

V 

volume  (m3) 

out 

system  outlet 

er 

electric  resistance  of  generation  of  heat  flux 

Greek  letters 

relat 

relative 

a 

seebeck  coefficient  (V/K) 

room 

ambient 

AV 

difference  of  electric  power  (V) 

negligible,  which  gives  an  approximation  closer  to  reality.  The 
model  inputs  are  only  the  room  temperature  and  the  heat  flux  sup¬ 
plied  by  the  residual  energy  source;  it  is  not  necessary  to  obtain 
the  parameters  experimentally. 

The  heat  exchangers  of  the  hot  and  cold  side  are  included  in  our 
model;  and  therefore,  the  temperatures  and  heat  fluxes  in  all  the 
components  of  the  thermoelectric  system  can  be  obtained  as  out¬ 
puts.  Another  advantage  is  that  our  model  determines  the  transi¬ 
tory  state  as  well  as  the  steady  state. 

2.  Objectives 

The  objectives  we  have  proposed  are: 

•  A  computational  model  design  capable  to  simulate  the  thermo¬ 
electric  generator  operation,  including  the  heat  transfer  from  the 
thermal  sources  to  the  faces  of  the  Peltier  module.  We  will  keep 
in  mind  the  Thompson  coefficient  and  all  the  thermoelectric 
properties  will  be  a  function  of  the  temperature. 

•  Validation  of  the  computational  model  with  experimental  data 
obtained  from  the  test  bench. 

•  Theoretical  and  experimental  study  of  the  influence  on  the  ther¬ 
moelectric  power  generated  of  the  following  parameters: 

-  Residual  heat  flux  and  electric  load  resistance  for  a  constant 
room  temperature. 

-  Room  temperature  and  electric  load  resistance  for  a  constant 
residual  heat  flux. 

-  Room  temperature  and  electric  load  resistance  for  a  constant 
temperature  gap  between  the  faces  of  the  Peltier  module. 

3.  Computational  model 

In  thermoelectric  devices,  for  both  refrigeration  and  generation 
applications,  the  temperatures  of  the  faces  of  the  Peltier  modules 
and  in  the  internal  nodes  can  not  be  determined  analytically  with¬ 
out  knowing  the  heat  flux  due  to  the  thermoelectric  effects.  The 


system  of  equations  in  order  to  calculate  the  temperatures  is 
nonlinear. 

The  model  solves  these  nonlinear  systems,  formed  by  the  ther¬ 
moelectricity  equations  and  the  heat  transfer  equations,  by  the  fi¬ 
nite  difference  method,  which  calculates  the  temperature  at 
different  points  separated  in  space  by  a  discrete  distance.  In  the 
transient  state,  the  temperatures  of  these  points  are  calculated  at 
discrete  periods  of  time.  With  this  purpose,  a  finite  period  of  time 
is  chosen  and  the  temperatures  for  all  the  points  are  recalculated  at 
the  end  of  this  time  interval. 

While  solving  the  model  using  the  implicit  finite  difference 
method,  the  values  of  heat  flux  can  be  determined  using  the  values 
of  the  temperatures  of  the  time  step  before. 

The  inputs  of  the  model  are:  the  geometric  data  and  the  material 
properties  of  the  elements  of  the  thermoelectric  system  studied 
(thermal  conductivity,  electric  resistivity,  Seebeck  coefficient  and 
specific  heat),  and  the  thermal  energy  value  supplied  to  the  sys¬ 
tem.  After  the  simulation,  the  outputs  of  the  model  are:  the  values 
of  efficiency,  electric  voltage  and  current,  electric  power  generated, 
temperatures  and  heat  flux  generated. 

Previous  works,  showed  in  Refs.  [3-13],  needed  the  tempera¬ 
ture  values  on  the  Peltier  module  faces  as  inputs  for  their  models. 
These  were  experimentally  measured  in  the  thermoelectric  gener¬ 
ation  system  for  each  heat  flux  supplied. 

Our  model  can  determine  all  the  heat  fluxes  and  temperatures 
of  the  thermoelectric  generator  components  with  only  the  heat 
flux  supplied  by  the  residual  energy  source  and  the  room  temper¬ 
ature  as  inputs.  This  model  does  not  need  previous  experimental 
values  of  the  device  temperatures. 

3.1.  Hypothesis  of  the  model 

The  materials  used  in  the  test  bench  shown  in  Fig.  2  are  consid¬ 
ered  isotropic. 

The  model  has  been  built  neglecting  the  magnetic  field  effects, 
so  the  Hall,  Nernst,  Ettingshausen  and  Righi-Ludec  effects  are  ne¬ 
glected  as  well. 
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3.2.  Equations  of  the  model 

The  model  solves  the  equation  of  heat  transfer  in  transitory 
state  for  one  dimension  and  one  dimensional  flow: 

pcpdT/dx  =  k(S2T/Sx2)  +  q*  (1) 

The  equation  of  heat  conduction  (1)  multiplied  by  the  volume  and 
applied  to  the  ith  node  as  a  function  of  the  thermal  resistances 
and  capacities  is: 

(Tit  ~  T'd/Rt-u  +  (T'j  -  T\)/Rij  +  Qi  =  (Tj  -  T,)C,/Sr  (2) 

The  generation  or  absorption  of  heat  related  to  node  i  is  represented 
by  Qit  and  its  expression  is  given  by  Eqs.  (6)-(9). 

The  thermal  resistance  between  nodes  i  and  j,  Eq.  (3),  and  the 
thermal  capacity  of  node  i,  Eq.  (4),  are: 

Rij  =  UjUkAi)  (3) 

Q  =  ViPicp  (4) 

Therefore,  the  thermal  resistance  and  capacity  of  the  Peltier 
module  is  studied  as  in  [18];  in  our  case,  we  divide  the  module  into 
ten  nodes,  and  we  get  a  thermal  resistance  between  nodes  and  a 
thermal  capacity  associated  with  each  node  as  a  function  of  its 
temperature. 


Qh  —  — N2  (XhlTh  +  I2Rcont  (6) 

Qx=N20CcITc  +  I2Rcont  (7) 

Qj  =  I2R0  =  NI22pL/A  (8) 

Qa  =  cr/AT  (9) 


The  temperature  of  the  semiconductor  terminals  is  used  to  calcu¬ 
late  the  heat  fluxes  due  to  the  Peltier  effect.  The  contribution  of 
the  contact  thermal  effects  was  developed  in  [2]. 

The  discretization  of  the  thermoelectric  generator  uses  symbols 
of  an  electric  analogy,  shown  in  Fig.  1.  The  model  assigns  different 
nodes  to  the  Peltier  module  in  the  ceramic  part  and  in  the  end  of 
the  semiconductor,  and  thus,  it  is  not  necessary  to  correct  the  tem¬ 
perature  gap  between  them.  This  is  an  improvement  to  the  analytic 
solution  of  [2],  which  is  necessary  when  only  the  experimental 
temperature  of  the  ceramic  part  of  the  Peltier  module  is  avail- 
able.With  the  scheme  shown  in  Fig.  1  and  Eq.  (5),  we  obtain  the  fol¬ 
lowing  system  of  equations: 

[M]  [T'j]  =  [T'j]  +  dz/Cj  [Qj]  (10) 

The  matrix  system  (10)  is  developed  as  a  function  of  the  ther¬ 
mal  resistances  and  capacities  of  the  device,  as  shown  in  Eq. 
(11): 
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(ii) 


For  the  contact  resistances  between  the  elements  of  the  test 
bench,  the  works  of  Ritzer  and  Lau  [21]  were  used. 

The  thermal  resistance  and  capacity  of  the  heat  extender  are 
calculated  using  Eqs.  (3)  and  (4),  giving  the  following  values: 

H-Heat  Extender  —  0.1677  [K/W] 

C Heat  Extender  —  211  [I/kg  K] 


In  our  case,  the  thermal  resistance  of  the  dissipater  and  the  insula¬ 
tor  related  to  the  ambient  are  calculated  experimentally  as  de¬ 
scribed  below: 

Grouping  the  equation  parameters  as  functions  of  the  node 
temperatures  and  the  heat  flux,  we  obtain  the  equation: 

(CjRj_ii)  +  T^T/CKl/Ki-i;  +  1  /Rij)  +  1]  -  T'jdz/CiRij 

=  Ti  +  QiSz/Ci  (5) 

The  model  incorporates  the  equations  from  the  thermoelectric  ef¬ 
fects.  These  effects  are  the  heat  flux  from  the  Peltier  effects,  Eqs.  (6) 
and  (7),  Joule  Eq.  (8)  and  Thomson  Eq.  (9),  with  no  magnetic  fields. 


The  Seebeck  coefficient  is  a  parameter  variable  with  the  tem¬ 
peratures,  so  if  it  is  included  in  the  model  as  a  function  of  the 
temperature,  we  get  the  following  expressions  for  the  electro¬ 
motive  force  generated  due  to  the  Seebeck  effect,  the  voltage 
in  the  ends  of  the  generator,  the  current  and  the  electric 
power. 


(oChTh  McTc)  ^  ^  I'i+i) 

i= 1 

AV  =  (m/1  +  m)Et 
I  =  £t/(l  +  rri)Ro 
Pout  =  a  VI 


(12) 

(13) 

(14) 

(15) 


4.  Experimental  work 

The  assembly  of  the  test  bench,  Fig.  2,  has  been  designed  for  this 
work,  and  it  is  composed  of: 
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Room 

Electrical 

Heat 

Ceramic 

Temperature 

Resistance 
&  Isolation 

Extender 

Layer 

Ceramic 

Cold 

Room 

Layer 

Dissipater 

Temperature 

_ ^ 

_ 

rrsr 


Sr.  St.  St_  °iot_ SiT.Sfr  0ijf°i4ri 


Fig.  1.  Scheme  of  the  thermal-electric  analogy  of  the  computational  model. 


Fig.  2.  Scheme  of  the  test  bench. 


•  A  calibrated  electric  resistance  that  simulates  the  thermal 
energy  source,  providing  a  heat  flux  Qer  with  a  maximum  value 
of  50  W. 

•  A  heat  extender,  with  a  known  thermal  resistance  that  commu¬ 
nicates  the  electric  resistance  with  the  Peltier  module. 

•  A  Peltier  module  Marlow  DT12-6L  type. 

•  A  dissipater  for  the  cold  side  of  the  Peltier  module,  with  a  ther¬ 
mal  resistance  experimentally  determined. 

•  A  decade  box,  CAM  METRIC  R420,  in  order  to  be  able  to  simulate 
with  different  load  resistances. 

The  principal  element  of  study  is  the  Peltier  module,  Marlow 
DT12-6L,  used  generally  in  thermoelectric  refrigeration.  We  have 
chosen  this  module  to  take  advantage  of  residual  thermal  energies 
of  low  temperature,  (lower  than  85  °C),  as  its  temperature  range  is 
the  most  suitable. 

In  order  to  determinate  with  accuracy  the  thermal  resistance  of 
the  dissipater,  we  made  some  tests  in  a  standardized  climatic 
room.  The  tests  were  made  with  a  prototype  that  used  a  known 
heat  flux  that  went  through  the  dissipater  with  the  same  temper¬ 
ature  room  conditions  as  those  of  the  thermoelectric  generator. 

In  Fig.  2,  it  can  be  seen  that  the  heat  flux  that  goes  through  the 
heat  extender  is  the  same  as  the  heat  absorbed  by  the  Peltier  mod¬ 
ule  (bn.  This  is  calculated  with  the  temperature  gap  between  the 
ends  of  the  heat  extender  and  the  thermal  resistance  of  the  heat 
extender  using  Eq.  (16). 

Q.in  —  AT Heat  Extender  / ^Heat  Extender  (16) 

In  an  open  circuit  case  (RL  =  0),  there  is  no  Peltier  effect,  and  thus, 
there  is  no  electric  power  generated  in  the  Peltier  module  (Pout  =  0), 


Table  1 

Thermal  resistance  of  the  dissipater. 


so  the  heat  flux  absorbed  is  Qin,  and  fits  with  the  flux  through  the 
dissipater  QDiSsiP;  With  this  data  and  experimentally  measuring 
the  temperature  gap  between  the  dissipater  and  the  ambient,  we 
determined  the  thermal  resistance  of  the  dissipater  using  Eq.  (17). 

R Dissip  —  AT Dissip- Ambient  /  Q.Dissip  (17) 

In  Table  1,  the  values  of  the  thermal  resistance  for  different  heat 
fluxes  of  the  dissipater  ClDiSsip  are  shown.  The  data  of  the  thermal 
resistance  of  the  dissipater  are  independent  of  the  thermoelectric 
device  installed,  and  the  thermal  resistance  of  the  dissipater  is  an 
input  data  for  the  model. 

The  temperature  probes  are  thermocouples  of  K  type,  and  the 
data  logging  is  ALMEMO  5590-2.  A  thermocouple  is  placed  in  each 
point  of  study  shown  in  Fig.  2.  The  tests  were  placed  in  a  climatic 
room,  CLIMATS  1440H  60/3,  to  keep  the  room  temperature 
constant. 

The  heat  flux  due  to  the  leaks,  Eq.  (19),  is  the  difference  between 
the  heat  fluxes  generated  by  the  electric  resistance,  Eq.  (18),  and 
that  absorbed  by  the  Peltier  module,  Eq.  (16). 

Qer  =  A  VerIer  (18) 

Qleaks  —  Qer  —  Qin  (19) 

The  thermal  resistance  between  the  electric  resistance  and  the 
ambient  is  determined  using  the  same  methodology  that  was  used 
to  calculate  the  thermal  resistance  of  the  dissipater,  using  two 
thermocouples,  one  placed  in  the  ambient  and  the  other  in  the 
electric  resistance.  Then  the  temperature  gap  between  the  electric 
resistance  and  the  ambient  is  known.  Thus,  with  the  value  of  the 
heat  flux  that  is  transmitted,  we  calculate  the  thermal  resistance 
as  is  shown  in  Table  2  for  different  values  of  the  incoming  flux. 

In  order  to  study  the  electric  power  generated,  it  is  necessary  to 
keep  in  mind  the  load  resistance  of  the  device  joint  to  the  thermo¬ 
electric  generator.  To  simulate  this  load  resistance  a  decade  box 
has  been  used,  CAM  METRIC  R420,  which  provides  a  variable  resis¬ 
tance  from  0.01  Cl  to  100  Cl. 
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Table  2 

Thermal  resistance  between  the  thermal  generator  source  and  the  room. 


5.  Results  and  discussion 


The  experimental  data  obtained  from  the  test  bench  are  com¬ 
pared  with  the  values  from  the  computational  model.  With  this 
comparison,  we  make  a  validation  of  the  model  to  determine  the 
error  between  the  experimental  and  the  simulated  values. 

The  sample  mean  Pout,e  of  the  measurement  results  is  estimated 
from  10  independent  observations  of  Ps  obtained  under  the  same 
conditions  of  measurement. 
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The  standard  uncertainty  u(P0Ut>£)  associated  with  Pout  e  is  the  esti¬ 
mated  standard  deviation  of  the  mean,  given  by  the  expression. 
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Therefore,  the  relative  uncertainty  of  the  experimental  mea¬ 
surement  is  estimated  as  follows: 


Ure,at(Pout,e)  =  (22) 

J  out, 8 

The  total  relative  uncertainty  of  the  experiments  was  acceptable.  It 
was  lower  than  1.5%  in  the  case  of  output  power  and  lower  than  2% 
in  the  case  of  temperatures. 

In  order  to  study  the  output  power  of  the  thermoelectric  gener¬ 
ator  system,  we  have  analyzed  three  different  cases: 


(I)  Constant  room  temperature,  varying  the  residual  heat  flux 
and  the  load  electric  resistance. 

(II)  Constant  residual  heat  flux,  varying  the  room  temperature 
and  the  load  electric  resistance. 

(III)  Constant  temperature  gap  between  the  Peltier  module,  vary¬ 
ing  the  room  temperature  and  the  load  electric  resistance. 


The  model  we  present  in  this  work  can  determine  the  temper¬ 
atures  and  the  heat  fluxes  of  all  nodes  of  the  thermoelectric  gener¬ 
ation  system  components.  It  only  needs  as  input  the  heat  flux 
supplied  from  the  hot  source  and  the  room  temperature.  It  is  not 
necessary  to  know  previously  the  electric  current  as  in  [16],  or 
[17].  In  our  case,  the  current  is  not  an  input  parameter  but  an  out¬ 
put  parameter  related  with  all  the  unknowns  by  Eq.  (14);  more¬ 
over,  it  uses  the  thermoelectric  parameters  as  a  function  of 
temperature. 

5.1.  Case  I:  Constant  room  temperature  varying  the  residual  heat  flux 
and  the  load  electric  resistance 


5 A  A.  Computational  model  validation 

In  the  tests  developed  for  the  present  case,  the  heat  flux  was 
constant.  This  heat  flux  simulates  the  residual  thermal  energy. 
The  electric  power  obtained  is  measured  for  each  of  the  load  resis¬ 
tances  tested. 
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Fig.  3.  Comparison  of  the  power  generated  from  the  model  and  the  experimental 
data  for  a  heat  flux  of  5  W. 
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Fig.  4.  Comparison  of  the  electric  power  generated  by  the  model  and  experimental 
data  for  a  supplied  heat  flux  of  30  W. 


Figs.  3  and  4  show  a  comparison  of  the  model  and  the  experi¬ 
mental  values  for  residual  heat  fluxes  of  Qin  =  5W  and  Qin  = 
30W,  respectively.  The  values  are  acceptable,  with  an  error  lower 
than  5%. 

The  maximum  electric  power  is  produced  for  values  of  load 
resistance  equal  to  the  internal  resistance  of  the  Peltier  module. 
This  behavior  was  demonstrated  analytically  in  [2].  For  a  constant 
room  temperature,  the  greater  is  the  supplied  residual  heat  flux, 
the  higher  is  the  temperature  of  the  Peltier  module. 

5 A  2.  Output  thermoelectric  power  -  experimental  and  simulated 
Our  model  has  into  account  the  variation  of  the  electric  resis¬ 
tance  of  the  Peltier  module  with  temperature.  Thus,  the  load 
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electric  resistance  that  gives  the  maximum  power  increases,  with 
the  temperature,  as  the  internal  electric  resistance  of  the  Peltier 
module  increases  with  the  temperature  as  well.  As  an  example,  it 
can  be  seen  that  if  a  heat  flux  of  30  W  is  supplied,  the  load 
resistance  value  for  the  maximum  power  generated  is  3.9  Q,  and 
if  the  heat  flux  has  a  value  of  5  W,  the  load  resistance  increases 
only  to  3.3  Q,  which  makes  the  rise  18%.  This  effect  is  shown  in 
Fig.  5. 

With  the  analytical  solution  of  the  ideal  thermoelectric  couple 
developed  from  Refs.  [3-1 1  ]  or  in  the  models  that  use  average  val¬ 
ues  of  the  thermoelectric  parameters,  Refs.  [12,13],  the  power  and 
efficiency  curves  are  calculated  as  a  function  of  the  electric  current 
generated,  keeping  constant  the  temperature  gap  between  the 
faces  of  the  Peltier  module. 

With  the  models  of  the  mentioned  works  it  is  not  possible  to 
determine  the  effect  shown  in  Fig.  5.  As  a  matter  of  fact,  increasing 
the  supplied  heat  flux  causes  the  absorbed  heat  due  to  the  Peltier 
effect  to  increase.  Thus,  the  temperature  gap  between  the  Peltier 
module’s  faces  increases.  In  the  cases  of  the  works  from  Refs.  [3- 
1 1  ],  it  would  be  necessary  to  calculate  experimentally  the  temper¬ 
ature  of  the  faces  of  the  Peltier  module  each  time  the  heat  flux  var¬ 
ies  in  order  to  determine  the  generated  thermoelectric  power. 


Fig.  5.  Thermoelectric  generated  power  maximum  for  each  curve  of  heat  flux 
constant  with  a  room  temperature  of  Troom  =  273  I<. 
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Fig.  6.  Temperatura  gap  between  the  faces  of  the  Peltier  module  as  a  function  of  the 
load  resistance  RL  for  a  constant  supplied  heat  flux. 


5.1.3.  Temperature  gap  between  the  sides  of  the  Peltier  module  - 
experimental  and  simulated 

The  model  calculates  the  temperatures  and  the  heat  fluxes  of  all 
the  elements  of  the  thermoelectric  device.  As  can  be  seen  in  Fig.  6, 
for  a  constant  heat  flux  supplied,  when  the  load  resistance  RL  in¬ 
creases,  the  temperature  gap  between  the  faces  of  the  Peltier  mod¬ 
ule  increases  (AT)  too.  The  simulation  results  with  the  model  have 
very  good  accuracy,  with  errors  lower  than  5%  with  the  experimen¬ 
tal  results.  As  an  example,  let  us  see  what  happens  for  a  supplied 
heat  flux  of  30  W  where  the  temperature  gap  between  the  faces 
of  the  Peltier  module  increases  10  K  (which  makes  an  increase  of 
27%)  when  the  load  resistance  is  varied  from  2  Q  to  16  Q. 

5.2.  Case  II:  Constant  waste  heat  flux  varying  the  room  temperature 
and  the  electric  load  resistance 

5.2.1.  Output  thermoelectric  power  experimental  and  simulated 

For  the  tests  made  for  a  constant  heat  flux  of  9.97  W  and  vary¬ 
ing  the  room  temperature  from  -25  °C  to  50  °C,  Fig.  7  shows  that 
the  experimental  values  obtained  of  generated  power  vary  with 
the  room  temperature.  This  fact  is  due  to  the  variation  of  the  ther¬ 
moelectric  properties  as  a  function  of  temperature. 
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Fig.  7.  Thermoelectric  generated  power  for  a  room  temperature  constant  with  a 
supplied  heat  flux  of  9.97  W. 
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Fig.  8.  Maximum  power  for  each  curve  of  constant  room  temperature,  with  a 
supplied  heat  flux  of  9.97  W. 
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Fig.  9.  Thermoelectric  generated  power  for  a  constant  room  temperature,  with  a 
temperature  gap  between  the  faces  of  the  module  of  20  °C. 


If  we  suppose  that  the  thermoelectric  properties  do  not  depend 
on  the  temperature,  for  the  tests  where  the  heat  flux  was  constant, 
the  value  of  the  electric  power  generated  would  be  constant,  inde¬ 
pendent  of  the  room  temperature.  Experimentally,  we  have 
checked  that  this  supposition  is  not  correct,  see  Fig.  7. 

The  model  and  experimental  values  are  shown  in  Fig.  7,  where 
it  can  be  checked  that  the  errors  are  lower  than  5%. 

Fig.  8  shows  the  maximum  power  curve  (for  a  constant  room 
temperature,  Fig.  7).  The  maximum  of  the  power  curve  is  produced 
for  a  load  resistance  that  is  greater  as  the  room  temperature  is 
greater.  For  this  case  of  study,  the  curve  of  maximum  power  has 
a  parabolic  shape,  reaching  the  highest  value  for  a  room  tempera¬ 
ture  of  -1  °C.  The  influence  of  the  temperature  on  the  thermoelec¬ 
tric  properties  depends  on  the  material  used,  so  with  other  Peltier 
modules,  different  curves  would  be  obtained. 

5.3.  Case  III:  Constant  temperature  gap  between  the  faces  of  the  Peltier 
module  varying  the  room  temperature  and  the  electric  load  resistance 

5.3 A.  Output  thermoelectric  power  experimental  and  simulated 

In  order  to  check  the  influence  of  the  room  temperature,  the 
cases  where  the  temperature  gap  between  the  faces  of  the  Peltier 
module  is  constant  have  been  simulated.  The  electric  power  gener¬ 
ated  is  shown  in  Fig.  9,  for  a  constant  temperature  gap  of  20  °C, 
with  varying  the  load  resistance. 

If  we  suppose  that  the  thermoelectric  properties  are  not  a  func¬ 
tion  of  the  temperature,  for  the  tests  where  the  temperature  gap  of 
the  faces  of  the  Peltier  module  is  constant,  the  maximum  power 
generated  (calculated  using  the  model  of  an  ideal  thermocouple 
or  the  models  that  suppose  average  values  for  the  thermoelectric 
properties)  would  be  constant,  independent  of  the  room  tempera¬ 
ture.  However,  the  simulations  and  the  experimental  data  show 
that  this  is  not  correct. 

For  the  case  of  our  study,  Fig.  9,  a  decrease  of  9%  of  the  maximum 
power  generated  was  obtained  between  the  lowest  room  tempera¬ 
ture,  Troom  =  -20  °C,  and  the  highest,  Troom  =  40  °C,  for  a  constant 
temperature  gap  between  the  faces  of  the  Peltier  module  of  20  °C. 

For  lower  values  of  load  resistance,  the  generated  power  is 
greater  as  the  room  temperature  is  lower.  For  higher  values  of 
room  temperature,  the  generated  power  is  greater  as  the  room 
temperature  is  greater,  as  shown  in  Fig.  9. 

Our  simulations  allow  choosing  the  load  resistance  more  suit¬ 
ably  as  a  function  of  the  ambient  (room)  temperature  of  operation. 


6.  Conclusions 

•  A  complete  computational  model  has  been  developed  that  is 
capable  to  simulate  the  thermoelectric  generation  of  the  Peltier 
module.  It  solves  the  equations  of  thermoelectricity  and  the  heat 
transfer  phenomenon.  The  thermoelectric  parameters  are 
defined  as  functions  of  temperature,  which  allows  us  to  keep 
in  mind  the  Thomson  effect. 

•  The  computational  model  has  been  validated  using  experimen¬ 
tal  data  of  a  test  bench  for  different  room  temperatures  and  dif¬ 
ferent  heat  fluxes.  The  errors  are  lower  than  5%. 

•  Our  computational  model  has  important  advantages  relative  to 
other  calculation  methods  from  the  literature  such  as: 

-  The  model  determines  the  thermoelectric  power  generated 
for  any  boundary  conditions  of  operation  (supplied  heat  flux 
to  the  Peltier  module  and  room  temperature).  It  is  not  neces¬ 
sary  to  obtain  parameters  experimentally. 

-  It  solves  the  system  of  equations  determining  the  tempera¬ 
tures  of  all  the  elements  in  the  thermoelectric  system  and 
the  heat  fluxes  involved.  It  determines  the  transitory  state  as 
well  as  the  steady  state.  This  makes  the  model  very  useful  as 
a  design  tool  of  thermoelectric  generation  systems. 

-  The  influences  of  room  temperature,  residual  heat  flux  and 
electric  load  resistance  on  the  behavior  of  the  thermoelectric 
generation  system  were  studied. 

-  It  was  experimentally  proved  and  demonstrated  that  the  elec¬ 
tric  load  resistance  that  gives  the  maximum  power  varies 
with  the  supplied  heat  flux.  Specifically,  in  our  case  of  study, 
it  raises  18%  when  the  electric  load  value  is  increased  from 
3.3  Q  (5  W)  to  3.9  Q  (30  W).  This  effect  can  be  observed  in 
the  computational  model  as  well. 

-  For  a  constant  supplied  heat  flux,  it  has  been  verified  that 
the  maximum  power  obtained  is  a  function  of  the  room 
temperature  and  the  load  resistance.  For  the  case  of  a  sup¬ 
plied  heat  flux  of  9.97  W,  the  curve  of  the  maximum 
power  as  a  function  of  the  room  temperature  has  a  para¬ 
bolic  shape,  reaching  the  maximum  at  the  temperature  of 
-1  °C. 

-  For  a  constant  temperature  gap  between  the  faces  of  the  Pel- 
tier  module,  the  influence  of  the  room  temperature  on  the 
maximum  power  was  studied.  In  our  case,  a  decrease  of  9% 
of  the  maximum  power  generated  was  obtained  between 
the  lowest  simulated  room  temperature,  -20  °C,  and  the 
highest  room  temperature,  40  °C. 
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